V ' ,J X , 4Sctat?school 

llioiiilsHEY, CALIFORBIA 939«-8002 





NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



ENERGY DISTRIBUTION OF CERENKOV 
RADIATION 
FOR 

FINITE FREQUENCY INTERVALS 
by 

Thomas M. Wilbur 



June 1987 



Thesis Advisor J. R. Neighbours 



Approved for public release; distribution is unlimited. 



T 2344 11 




SECuRit’’ Cl ASS.FiCATlON OF TwTT PAGE 



REPORT DOCUMENTATION PAGE 



la REPOST SECURITY Classification 

UNCLASSIFIED 



ib RESTRICTIVE markings 



U SECURITY Classification auTmqRiTy 



ED DEClASSiFiCAT.ON i OOWNGRAOiNG SChEOULE 



3 OiSTRiBUTiON/ AVAILABILITY of report 

Approved for public release; distribution 
is unlimited 



4 PERFORM. NO ORGAN. EATlON REPORT NUM8E R(S) 



S MQNiTORiNG ORGANiEAT.ON REPORT NuV8£R(S) 



fee NAME OF PERFORMING ORGANISATION 

Naval Postgraduate School 



6d off.ce symbol 
(it 4ppn<*b>*) 

61 



>4 NAME OF MONITORING ORGANiSATiON 

Naval Postgraduate School 



cx aOORESS iGry Sure 4nd hP Code) 

Monterey, California 93943-5000 



7b aOORESS (Cry sure 4nd JiPCode) 

Monterey, California 93943-5000 



NAME OF FuNOiNG / SPONSORING 
organisation 



8b OFFICE SYMBOL 

(If 4ppliC4bl*) 



9 PROCUREMENT INSTRUMENT iQE N TiF iC A T ION NUMBER 



Be AOORESSiCiry Sure end S/P Cod*) 



10 SOURCE OF FuNOiNG NUMBERS 



PPOGAAM 


P»0,ECT 


TASK 


AORK jNir 


ELEMENT NO 


NO 


NO 


ACCESS ON NO 



i TlE (include SecL/r/fy CUwti(4tion) 

ENERGY DISTRIBUTION OF CERENKOV RADIATION FOR FINITE FREQUENCY INTERVALS 



PERSONA^ auThOR(S) 
Wilbur, Thomas M. 



)j OF REPORT 


l )D T ME COVERED 


14 DATE OF REPORT (Yter Month 04y) 


’S PAGE COuNT 


Master's Hies is 


FROM TO 


1987 June 


50 



6 supplementary notation 



COSATi COOES 


*€u-3— 


GROUP 


SuB-GPOUP 








i 


1 





18 SUBJECT TERMS (Continue on reverie it neceisery 4r>d identity O y block number) 

Cerenkov radiation 
Electron Accelerator 



1 aBs t RaCT (Continue on reverie it neceuety 4na • oennty by block number) 

The equation defining the energy radiated per unit solid angle due 
to Cerenkov radiation is analyzed in detail, including the effects of 
all equation variables for a hypothetical electron accelerator 
experiment. Specifically, various finite frequency intervals are used 
in an effort to determine the optimum means of determining the details 
of a charge bunch in a high energy electron accelerator. In particular, 
it is shown how narrowband measurements as a function of angle may yield 
information on both the beam path length and the bunch charge 
parameters. As an aid to the analysis, an interactive Fortran program 
is presented that allows for any specific experimental parameters, with 
options for various output types as desired. 



0 0 S '« BuTiON I availability OF ABSTRACT 
G) -NCLaSSiFiEDUNL'MiTEO □ SAME AS RPT O OTiC USERS 



21 abstract security classification 

UNCLASSIFIED 



2)4 NAME OF RESPONSIBLE iNDiViOUAI 



22b TELEPHONE (Include Artj Cobf) 

(408) 646-2291 



22 c OFFICE SYMBOL 

6 1 Nb 



00 FORM 1473. 84 mar 



83 apr edition be u»ed until circuited 
All other editions 4<e obsolete 



SECURITY CLASSIFICATION OF ThiS PAGE 



1 



Approved for public release; distribution is unlimited. 



Energy Distribution of Cerenkov Radiation 

for 

Finite Frequency Intervals 



by 



Thomas M. Wilbur 
Lieutenant, United States Navy 
B.S., Pennsylvania State University, 1978. 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN PHYSICS 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1987 



ABSTRACT 



The equation defining the energy radiated per unit solid angle due to Cerenkov 
radiation is analyzed in detail, including the effects of all equation variables for a 
hypothetical electron accelerator experiment. Specifically, various finite frequency 
intervals are used in an effort to determine the optimum means of determining the 
details of a charge bunch in a high energy electron accelerator. In particular, it is 
shown how narrowband measurements as a function of angle may yield information on 
both the beam path length and the bunch charge parameters. As an aid to the 
analysis, an interactive Fortran program is presented that allows for any specific 
experimental parameters, with options for various output types as desired. 
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THESIS DISCLAIMER 



y,&S0l 

a - 1 



The reader is cautioned that computer programs developed in this research may 
not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of computational and 
logic errors, they cannot be considered validated. Any application of these programs 
without additional verification is at the risk of the user. 
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I. INTRODUCTION 



A. HISTORY 

The existence of the phenomenon now known as Cerenkov radiation was 
observed as early as 1910, most notably by Madam Curie. However, during the 
ensuing years, there was other work being completed by those familiar with Madam 
Curie's observation that masked a detailed study of the phenomenon. Jelley [Ref. 1] 
explains their reasons in greater detail. In 1926, Mallet made the first deliberate 
attempt to study and explain “the very faint emission of a bluish-white light from 
transparent substances”. Experimentation concerning this phenomenon continued 
through the 1920's and 1930's by Cerenkov and Mallet, but a viable theory explaining 
the process would not be proposed until 1937 by Frank and Tamm. Their theory, for 
which they were consequently awarded the 1958 Nobel Prize in Physics, was found to 
be in excellent agreement with the experimental results obtained by Cerenkov. The 
advent of more sensitive light detectors and other experimental equipment accounted 
for more in depth studies of Cerenkov radiation in the 1940's and 1950's. [Ref. 1: pp. 
1 - 8 ] 

Since then, numerous experiments have been conducted and many papers and 
academic theses published that have helped to more clearly understand and explain the 
Cerenkov radiation phenomenon. Some of the more recent investigations include; 
periodic electron bunches of finite emission lengths, time developement of bunch 
charges and Cerenkov radiation in the x-ray region, (see Bibliography). 

B. BACKGROUND 

Work on this thesis was motivated by two factors. First, using previously derived 
equations and expressions describing the energy radiated per unit solid angle due to 
Cerenkov radiation, [Ref. 2] a fortran program was written to run on the IBM 3033 
mainframe. The program, (see Appendix), was written interactively to allow all 
parameters to be changed in order to tailor the output to fit any specific experimental 
setup. 

Second, using the aforementioned Fortran program, a series of outputs were 
generated to be used in comparing the theoretical results of previously derived 
expressions with those of actual experimentation. Since there are literally an infinite 
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number of combinations of various parameters, only a few were chosen to be included 
in this work. The parameters chosen are consistent with those expected while 
conducting experiments on any high energy electron accelerator. 
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II. THEORY 



A. CERENKOV EFFECT 

Cerenkov radiation can most easily be explained by describing the interaction of 
a single electron within a transparent medium. Consider a non-relativistic electron. 
While traversing the medium, the electron will tend to locally polarize the adjacent 
atoms in the medium, instantaneously creating a temporary dipole within the material. 
Therefore, as the electron moves through the medium, an electromagnetic pulse is 
generated. However, since there is complete symmetry’, there will be no net field 
generated. See f igure 2.1. 



Now suppose an electron is travelling at a speed comparable to or exceeding the 
speed of light for the medium. Although symmetry is preserved in the azimuthal plane, 
along the axis of motion a resultant dipole field is generated. These fields are set up at 
each element along the electron's track, radiating a brief electromagnetic pulse. 
Provided the electron is moving at a speed greater than the speed of light within the 
medium, the wavelets formed from all elemental positions on the track can be in phase 
and thus produce a resultant electromagnetic field. [Ref. 1: pp 3-6J 

Figure 2.2 depicts the relationship of the distance {BZ) travelled by an electron in 
time At versus the distance {BA) covered by an emitted electromagnetic pulse during 
the same time. The distance {BZ) travelled by the electron is given by: 




(A) 



(B) 



Figure 2.1 Polarized Atoms in a Dielectric. 
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1)7 = fk () At, 



(cqn 2.1) 



where p is the ratio of the electron speed to the speed of light in a vacuum ,c Q . The 
distance ( BA ) travelled by the electromagnetic wave is given by: 

BA = cAt, (cqn 2.2) 

where c is the speed of light for the medium. For a medium with an index of refraction 
n, c = c 0 /n. From equations 2.1 and 2.2, the “cercnkov relation*’ is obtained. 
| Ref. I: pp 6*S| 



cos0 c = 1/pn. 



(eqn 2.3) 



VAt 




Figqre 2.2 Illustration of Cerenkov Radiation. 

Due to the finite path length at which the radiation is detected, the Cerenkov 
cone angle is normally shifted away from 0 c . As the frequency at which the emitted 
radiation is measured or the distance from the source is increased, the radiation cone 
angle approaches 0 c . As this occurs, an increasing fraction of the total Cerenkov 
radiation is found at 0 c . [Ref. 2] 
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B. MATHEMATICAL INTERPRETATION 



1. Energy Equation 

As previously discussed, for a finite path length, diffracted Cerenkov radiation 
effects will be seen at angles other than the Cerenkov angle. Previous work on this 
phenomenon (Ref 2] has resulted in an expression for the emitted Cerenkov radiation. 
The energy radiated per unit solid angle within the frequency range dv by a single 
bunch charge q travelling a finite distance L is: 

E(v,£)dv = QR 2 dv, (eqn 2.4) 

where Q is a constant whose magnitude is dependent on the index of refraction and the 
total charge in Coulombs. 



Q = ncq 2 /8rc 2 . (eqn 2.5) 

The radiation function R is given by: 

R = 2;ujsin0I(u)F(X), (eqn 2.6) 

where 0 is the angle measured from the line of travel of the bunch charge to the 
direction of propagation of the emitted radiation, I(u) is the diffraction function, and 
F{k) is a form factor. For convenience, the length of travel of the bunch charge is 
measured in units of the wavelength of the emitted radiation within the medium. The 
dimensionless beam length parameter r\ serves this purpose and is defined as the ratio 
of the length L from the source to the wavelength X of the emitted radiation. 

x\ = L/X. (eqn 2.7) 

The diffraction function is defined as: 

I(u) = sin(u)/u, (eqn 2.8) 

where u is dependent on both the angle 0 and the beam length parameter r\. 

u = ;tn{(l/np) - cos0}. (eqn 2.9) 
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Equations 2.5 through 2.9 can then be used to redefine equation 2.4 in terms of 
experimental constants, (i.e. n, L, P etc.) and the single remaining variable of 
frequency. [Ref. 2] 

E(v,/c)dv = QG 2 F 2 (A)sin 2 (Av}dv. (eqn 2.10) 

For a given experimental setup, A and G are functions of the angle only and are 
defined as follows: 



A= {Lrc/c}»{( l/nP)-cos0}, and (eqn 2.11) 

G = {2sin0}/{(l/np)-cos0}. (eqn 2.12) 

The radiation patterns arising from equation 2.4 can be thought of as an 
oscillating sin 2 x function modulated by an envelope. If F(/c) is neglected, then G 2 acts 
as the modulating envelope. The actual form of the envelope will vary with nP, but 
will be constant for a given experimental setup. Figure 2.3 shows a plot of G 2 as an 
illustrative example for values of nP two percent above and below threshold, 
(np= 1.0 ±0.02). 

Finally, the dimensionless form factor F(/c) is related to the Fourier transform 
of the bunch charge distribution. The Fourier components of any bunch charge are 
defined as: 



p(£) = Jp(r)e^d 3 r, (eqn 2.13) 

where p(r) defines the charge distribution of the bunch charge. Once the Fourier 
components are determined from equation 2.13, the form factor is found by the 
relation: 



p(A) = qF(A). (eqn 2.14) 

Here we consider only a line charge with distribution p(z) which is travelling in the z 
direction. Thus, r-+z in equation 2.13 and &-+k z in equation 2.14. In this work, only 
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G (0) 



CERENKOV RADIATION ENVELOPE 




Figure 2.3 Illustrative Example of the Cerenkov Radiation Envelope. 
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the trapezoidal bunch charge distribution is used although many others are available 
and easily used within the Fortran program. [Ref. 3] 

2. Trapezoidal Form Factor 

The first step in obtaining any type of form factor is to determine the assumed 
geometry’ of the emitted bunch charge. Figure 2.4 depicts the geometry for the 
trapezoidal bunch charge. The functions defining the positive half of the trapezoid are 
given by: 



fj(x) = A, and 



(eqn 2.15) 



f 2 (x) = (A/(d-b)}*(d-x), (eqn 2.16) 

where fj defines the trapezoid from 0 to b and f 2 defines the trapezoid from b to d. 
The amplitude of the bunch charge is denoted as A with units of coulombs per meter. 
The values b and d have units of meters and are easily obtained from the pulse length. 




Figure 2.4 Trapezoidal Bunch Charge. 

The fourier components of this trapezoid are defined by equation 2.13 where 
p(k z ) is the sum of the two functions fj(x) andf 2 (x), (ie. the charge distribution of a 
single bunch charge), and k z is the component of the wave vector k. The details of 
solving equation 2.13, though tedious, are rather straightforward. [Ref. 3] 
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Equation 2.17 is the solution to equation 2.13 and provides the Fourier components 
in the case of the trapezoid. 

p(k z ) = { (2A)(d-b)/k 2 z }•{ cos(k z b) - cos(k z d) }. (eqn 2.17) 

The value of A is found by equating q, the total bunch charge, to the integral 
of the charge distribution. 

q = -oof° P( z ) dz = A(b + d). (eqn 2.18) 

By substituting the results of equations 2.17 and 2.18 into equation 2.14. the form 
factor for the trapezoidal charge distribution is: 

F(k z ) = { 2q/(d 2 -b 2 ) } • ( 1/k 2 } • ( cos(k z b)-cos(k z d) }. (eqn 2.19) 

If the medium is assumed to be non-dispersive, then co = ck z and remembering 
that k z = A'cosG, the form factor can be written in terms of the angular frequency 
2nv = co. After substituting, the explicit relation for the trapezoidal form factor is 
given by equation 2.20 as follows: 

F(k z ) = { c 2 /2;t 2 (d 2 -b 2 ) }•{ cos(Bv)-cos(Dv) )•{ l/v 2 }, (eqn 2.20) 
where B and D are given by equations 2.21 and 2.22 respectively. 

B = 2rcbcos0/c, and (eqn 2.21) 



D = 27tdcos0/c. (eqn 2.22) 

It is clearly seen that once the parameters for the bunch charge are defined, 
the form factor is highly dependent on the inverse square of the frequency v. The 
expression for the energy radiated per unit solid angle within the frequency range dv 
written in terms of frequency will then be equation 2.4 with the form factor as given in 
equation 2.20. [Ref. 4] 
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(eqn 2.23) 



Edv = [ QG 2 sin 2 (Av) ] 

•[ {c 2 /2;i 2 (d 2 -b 2 )}{cos(Bv)-cos(Dv)} { 1/v 2 } ] 2 dv. 

The radiated energy is a function of frequency through both the sin 2 (Av) term and the 
form factor. 
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III. RESULTS AND ANALYSIS 



The Fortran program written for this work has various output capabilities. The 
first type of output contains data equating the explicit value of E (eqn.2.4), in joule- 
seconds, to a specific frequency at a particular angle 0. Another output type provides 
data for graphing the radiation envelope (eqn.2.3) or the form factor (eqn. 2.20) used in 
this thesis. These graphs are useful in understanding the complex nature of the 
variables involved in solving equation 2.4. 

The most useful output type is that which solves, (i.e. integrates), equation 2.4 
over a predetermined frequency range for all desired angles from the beam path. A 
careful study of these graphical outputs will provide insight into the particulars 
necessary for any experimental setup so that the most accurate and informative results 
can be obtained. 

As previously discussed, the output obtained during this work are based on a set 
of parameters that uniquely define a specific electron accelerator experimental setup. 
The following input parameters (with their Fortran variable names), are required: 

1. Accelerator Beam Energy (EBEAM). 

2. Total Periodic Electron Bunch Charge (CUE). 

3. Index of Refraction (IND). 

4. Distance from Source to Detector (LENGTH). 

5. Cerenkov Radiation Frequency Interval (MIDNU.ENU). 

6. Incremental Change in Frequency Interval (DNU). 

7. Angle of Interest, as Measured from the Charge Path (DTHETA). 

8. Bunch Charge Form Factor (FORFAC). 

9. Trapezoidal Bunch Parameters (BEE, DEE). 

The specific values for each variable used to obtain the outputs in this work are 
as follows: 

1. Beam Energy - 50 MEV. 

2. Bunch Charge - .001 Coulombs. 

3. Index of Refraction - 1.000268 (air). 

4. Source to Detector Distance - 100 Meters. 

5. Frequency Interval - As noted. 

6. Frequency Increment - 0.5 MHz. 
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7. Solid Angle * As noted. 

8. Form Factor - Unity and Trapezoid. 

9. Trapezoidal Bunch Parameters. 

a. Top-50 nanoseconds. 

b. Base-60 nanoseconds. 

For the given index of refraction and beam energy, the “Cerenkov Relation” (eqn.2.3) 
gives the critical Cerenkov angle as 1.9°. 

Three basic assumptions were made in the course of this work. First, it is 
assumed that the shape of the trapezoidal bunch charge remains unchanged as it 
travels through the medium. This is reasonable since the distances involved are 
relatively short with respect to the speed at which the bunch charge is travelling. 
Second, the permeability of free space [i is used vice that of the actual medium of air. 

The final assumption is that the speed of the electron bunch remains unchanged 
throughout it's travel. Although there is some Bremstrahlung radiation emitted, it has 
a negligeable effect on the speed of the bunch charge and therefore P is assumed to 
remain constant. 

Figure 3.1 depicts the energy radiated per unit frequency at angles of 30, 45 and 
60 degrees using a form factor of unity. As would be expected from equation 2.23, a 
snrx function results since Q (eqn.2.5), A(eqn.2.11) and G(eqn.2.12) are all constant at 
the given angle 0, the only variable being the frequency v. The difference in 
amplitudes is due to the variation of G with 0. The difference in the periods is due to 
the function A, which, as the argument of the sin term, effects the periodicity. 

Once a form factor other than unity is used, the output pattern is changed 
significantly. Equation 2.23 shows the frequency dependence of the trapezoidal form 
factor. Figure 3.2 is the form factor obtained at 0 degrees over the frequency range 
10-100 MHz. The damping of the form factor waveform is due to the inverse square 
relationship of F(&) with frequency. Figures 3.3 and 3.4 depict the square of the form 
factor as it is used in equation 2.4. Different frequency intervals are used but each 
waveform is computed at the same angle, 0 = 0°. As with Figure 3.2, Figures 3.3 and 
3.4 each dampen with increasing frequency. The increased periodicity of Figure 3.4 is 
due primarily to the w’ider frequency range. At angles other than 0 = 0° the form 
factor will be similar in shape except scaled by changes in B and D due to the cos0 
factor. [Ref. 4] 
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ENERGY DISTRIBUTION (F(K)=1.0) 




Figure 3.1 Energy Radiated per Unit Frequency, (9 = 30°, 45° and 60°). 
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TRAPEZOIDAL FORM FACTOR 




Figure 3.2 Trapezoidal Form Factor, (0 = 0°). 
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TRAPEZOIDAL FORM FACTOR 




Figure 3.3 Square of the Trapezoidal Form Factor, 10-100MlIz., (0 = 0°). 
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. 01 * 



TRAPEZOIDAL FORM FACTOR 




Figure 3.4 Square of the Trapezoidal Form Factor, 200-800 MI lz., (0 = 0°). 
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Figure 3.5 is a plot similar to that of Figure 3.1 except that the unity form factor 
has been replaced with that of the trapezoid. As expected, the use of a specific form 
factor significantly effects the output. The output of interest will be the integration of 
all outputs similar to Figure 3.5 at each angle of interest. Before analyzing the 
integrated output with the trapezoidal form factor, it is useful to look at Figures 3.6 
through 3.8 w'hich depict the integration of equation 2.23 while maintaining unity form 
factor. As noted earlier, an energy peak at the Cerenkov angle of 1.9° is expected. 
Although the peak appears to be centered at about 10° in Figure 3.6, as the frequency 
interval increases from 10-100 MHz. to 10-10,000 MHz. the same curve shifts closer to 
the Cerenkov angle. Concurrently, the graphs become smoother as more data points 
are used at each interval 0. The aberration found in Figure 3.8 may be an artifact of 
the Fortran program. The extraordinarily large magnitude of Figure 3.8 is due to the 
existence of a pole in variable G (eqn.2.12) of equation 2.10 as 0 approaches the 
critical Cerenkov angle. However, since the variable A is zero at the Cerenkov Angle, 
equation 2.10 remains finite. [Ref. 2] 

Figure -3.9 shows the output of equation 2.23 integrated over a relatively 
broadband frequency range of 100 to 1,000 MHz. Figure 3.10 is the identical output 
plotted over a narrower range of 0 to enhance the output detail. In each case, the 
detail of the output is fuzzy at best, making it difficult to conduct any meaningful 
analysis. In an attempt to force more detail in the output, the frequency range over 
w'hich equation 2.23 is integrated was narrowed. Figure 3.11 depicts the energy spectra 
integrated over 20 MHz. centered around 100 MHz. All other parameters are equal to 
those used in Figures 3.9 and 3.10. Although similar in general form , Figure 3.11 
shows more detail. Then, changing to a narrower frequency band of 2 MHz. (again 
centered around 100 MHz.), Figure 3.12 depicts a clearly defined output that is similar 
to a frequency modulated carrier w r ave. Figure 3.13 show's the same output for a 
narrower range of 0 to show more detail. 

The Cerenkov radiation envelope discussed earlier( Figure 2.3) has an impact on 
the outputs of Figures 3.12 and 3.13, but once n0 is selected the envelope will effect 
only the amplitude of this and similar graphs. The null points of Figures 3.12 and 3.13 
can then be calculated by analyzing the components of the radiation function(eqn.2.6). 
The expression for the trapezoidal bunch charge form factor(eqn.2.20) can be expressed 
as a product of a constant and two sin functions using the following trigonometric 
identity. 
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E (X10 JOULE-SECONDS) 



ENERGY DISTRIBUTION/! (F(K)= TRAPEZOID) 




Figure 3.5 Energy Radiated per Unit Frequency, (0 = 30°, 45° and60°). 
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ENERGY (JOULES) 



CERENKOV(IO-IOOMHZ) 




Figure 3.6 Energy per Unit Solid Angle, KMOOMHz., (Unity Form Factor). 
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ENERGY (X10*JOULES) 



CERENKOV(IO-IOOOMHZ) 




l igure 3.7 Energy per Unit Solid Angle, lO-lOOOMHz., (Unity Form Factor). 
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ENERGY (XIO'JOULES) 



CERENKOV(IO-IOOOOMHZ) 




Figure 3.8 Energy per Unit Solid Angle, 10-10, 000MHz., (Unity Form Factor). 
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cosct - cosy = -2sin{ (a/2) + (y/2) }*sin{ (a/2)-(y/2) }. 



(eqn 3.1) 



After a great deal of algebraic manipulation, it can be shown that the form factor will 
be of the form: 

F(&) = const*. sin(£)*.sin(q>), (eqn 3.2) 

where £ = 7tcos0(d + b)v/c, and <p = 7icos0(d-b)v/c. 

For the trapezoidal bunch charge parameters used in this work, c/(d-b) equates 
to a value of 200 MHz. and c/(d+b) to a value of 18.2 MHz. Then, using the mean 
frequency of 100 MHz. found in Figure 3.13, one need only to equate each of the sin 
function arguments to integer values of n to mathematically determine the values of 0 
at which equation 3.2 w r ill equal zero. Table 1 lists the values of 0 at which the form 
factor will force zeroes in the energy equation. 



TABLE 1 


VALUES OF 0 AT WHICH THE FORM FACTOR EQUALS ZERO. 


sin£ = 0 


sintp = 0 


0j = 79.5°. 


undefined. 


0 2 = 68.7°. 


undefined. 


0 3 = 56.9°. 


undefined. 


0 4 = 43.3°. 


undefined. 


0 5 = 24.5°. 


undefined. 



At the mean frequency of 100 MHz., sin (p = 0 cannot be determined since it 
requires solving for an angle 0 w r hose cos is greater than unity. However, by 
comparing the values of 0 in column 1 of Table 1 to Figure 3.13, one sees a correlation 
of the null points to what can be considered the envelope of the output. It is clear 
then that the form factor acts as the modulator of the energy spectra of the Cerenkov 
radiation. 
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The "carrier wave" nulls can be determined in a similar manner by analyzing the 
diffraction function (eqn.2.8) found in equation 2.6.. Equating u to integer values of k 
and solving for 0 will provide all angles 0 at which the energy output will be zero. 
Table 2 lists a sample of the results for which sinu equals zero. Again, comparing the 
values of 0 in Table 2 to Figure 3.13 , one sees that the calculated nulls equate to those 
of the "carrier wave". 



TABLE 2 

VALUES OF 0 AT WHICH SINU EQUALS ZERO. 

sinu = 0. 

0 8 = 40.6°. 

0 9 = 43.1°. 

0 10 = 45.6°. 

0jj = 47.9°. 

0 12 = 5O.2°. 

0 13 =52.4 o . 

0 h =54.6°. 

0 J5 = 56.6°. 

0 I6 =58.7*. 

0 17 = 6O.6 o . 



Figures 3.14 and 3.15 are graphs of the Cerenkov energy output over the same 
narrowband frequency interval of 2 MHz., but at the higher frequency range of 499 to 
501 MHz. Although some detail is lost compared to the outputs of Figures 3.12 and 
3.13, the same analysis can be conducted to determine the null points of both the 
"carrier wave" and the modulating form factor. Tables 3 and 4 list the values of 0 for 
each null point similar to those found for Figures 3.12 and 3.13. The range of 0 was 
chosen so that a meaningful comparison could be made with the outputs of Figures 
3.14 and 3.15. 
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ENERGY (X10 3 JOULES) 



CERENKOV (100-1000MHZ) 




Rgure 3.9 Energy per Unit Solid Angle, lOO-lOOOMIIz., (Trapezoidal Form Factor). 
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ENERGY (JOULES) 



CERENKOV (100-1000MHZ) 




Figure 3.10 Energy per Unit Solid Angle, 100-1000MHz., (Trapezoidal Form Factor). 
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ENERGY (JOULES) 



CERENKOV (90-110MHZ) 




Figure 3.11 Energy per Unit Solid Angle, 90-110MlIz., (Trapezoidal Form Factor). 
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ENERGY (JOULES) 



CERENKOV (99-101MHZ) 




Figure 3.12 Energy per Unit Solid Angle, 99- 101 MHz., (Trapezoidal Form Factor). 
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ENERGY (JOULES) 



CERENKOV (99-101MHZ) 




Figure 3.13 Energy per Unit Solid Angle, 99-101MHz., (Trapezoidal Form Factor). 
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TABLE 3 

VALUES OF 0 AT WHICH THE FORM FACTOR EQUALS ZERO. 



sin£ = 0. 


sin(p = 0. 


0j = 87.9°. 


0j = 66.4°. 


0 2 = 85.8°. 


0 2 = 36.9°. 


0 3 = 83.7°. 


undefined. 


0 4 = 81.6°. 


undefined. 


0 5 = 79.5°. 


undefined. 


0 6 = 77.4°. 


undefined. 


0 ? = 75.2°. 


undefined. 



TABLE 4 

VALUES OF 0 AT WHICH SINU EQUALS ZERO. 



sinu = 0. 



e 

e 

e 

0 

o 

o 

0 

0 



148 

149 

150 

151 

152 

153 

154 

155 



83.6°. 

83.9°. 

84.3°. 

84.6°. 

84.9°. 

85.3°. 

85.6°. 

85.9°. 
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ENERGY (JOULES) 



CERENKOV (499-501MHZ) 




Figure 3.14 Energy per Unit Solid Angle, 499-501 MI lz., (Trapezoidal Form Factor). 
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ENERGY (JOULES) 



CERENKOV (499-501MHZ) 




Figure 3.15 Energy per Unit Solid Angle, 499-501 MHz., (Trapezoidal Form Factor). 



IV. DISCUSSION 



Figures 3.6 through 3.8 show that the energy radiated per unit solid angle from a 
charge distribution with unity form factor tends to be concentrated into an ever 
narrowing direction as the range of integration over frequency is increased. In the limit 
of all frequencies, the radiation is confined to the Cerenkov cone. The effect of 
including all frequencies is the same as allowing the beam path length to go to infinity, 
and the resulting radiation patterns are characteristic of that from a point charge. 

The zeroes in the radiation patterns for trapezoidal charge distributions, Figures 
3.9 through 3.15, are tabulated in Tables 1 through 4. Since the tables and figures are 
calculated from the same theory, they are, of course, in agreement. Tables 2 and 4 
show the zeroes for the "carrier" oscillation which arises from the length of the electron 
beam path from source to detector. Tables 1 and 3 show the zeroes in the c and (p 
oscillations resulting from the length (d+b) and rise time (d-b) respectively of an 
individual charge bunch. Both Tables show zeroes dictated by the length of the bunch 
charge. Only two zeroes from the pulse rise time occur; both at the higher frequency 
of measurement. 

In the cases studied here, the bunch length (d + b) equals 16.5 meters and the 
pulse rise (d-b) equals 1.5 meters. At the measuring frequency of 100 MHz. the ratio 
of bunch length to wavelength is 5.5; large enough to resolve the details of the bunch 
size. Even at the higher frequency of 500 MHz. the ratio of pulse rise to wavelength is 
2.5 and consequently the rise and fall of the pulse have little effect on the radiated 
output. [Ref. 5] 
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V. CONCLUSIONS AND RECOMMENDATIONS 



After studying the figures depicting the radiated Cerenkov energy in this work, 
one could initially deduce that the use of narrowband detectors would be optimum for 
detecting radiation during any electron accelerator experiment. As shown and 
previously discussed, the information produced is more detailed when narrowband 
frequency ranges are used compared to wider band measurements. In particular, if one 
frequency band is considered as a function of angle, the calculated radiation energy has 
much detail as displayed in Figures 3.12 and 3.13 The structure depends on both the 
path length and bunch length, and, in principle, may yield information about these 
parameters. One must also consider the logistics of any experimental setup and any 
engineering problems involved to yield the maximum information. 

For the broadband output with a unity formfactor, the energy output has a peak 
at the critical Cerenkov angle as expected. This output is typical of a point charge. 
However, since the nominal electron accelerator bunch charge output has a finite 
distribution, there is little information to be extracted from this particular output. 

When the effect of the form factor is included, there exists a definite Cerenkov 
radiation peak perpendicular to the path the electron bunch is travelling. This 
phenomenon may be significant with respect to the optimum position at which the 
maximum energy produced by the radiation can be detected. 

If a narrowband detection system is realistically feasible, it can be used to 
determine all nulls that arise in the energy equation due to the diffraction function and 
the form factor. Although narrowband detection will produce the most detailed 
output, its usefulness may be limited by the availability of equipment designed to 
discriminate the output within the desired frequency band limits. 

In order to obtain significant information about the details of the charge 
distribution within a bunch, it will be necessary to perform experiments at frequencies 
such that the ratio of rise time to frequency is at least 5. For real beams, where the 
rise times are of the order of a few nanoseconds, this means that narrowband 
measurements must be performed at microwave frequencies. 
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To augment this research the following is recommended: 

1. Analyze outputs with form factors other than trapezoidal, (i.e. 
triangular etc.). 

2. Determine the optimum (hieher) frequencv bands for various 
intervals at which useful information can be extracted. 



rounded, 

frequency 
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APPENDIX 
FORTRAN PROGRAM 



C PROGRAM TO COMPUTE THE ENERGY RADIATED PER UNIT SOLID ANGLE TRI00010 

C OVER A SPECIFIC FREQUENCY RANGE 
C 

C * ERG ARRAY SIZE MUST CORRESPOND TO THE NUMBER OF 

C * INCREMENTS OF THETA 

C 

DIMENSION ERG( 1800 ) 

C 
C 

REAL ERG, COSANG, COSARG, COSGAR, KKONST , DEE ,BEE 
10 REAL FORFAC,EBEAM, LENGTH, KONST, C, GEE, TOP, BOT 
20 REAL ATHETA,RADFUN, PI, EREST, GAMMA, THETA, IND,MIDNU 

30 REAL CONST, MU, CUE, TERG ,TTERG,DNU,BNU,DEGINC ,DATEND 

INTEGER IF0RM,INCREM,DL00P,D2L00P, SELECT, NUMBER 
CHARACTER*! RESPON,ANSW 
C 
C 

C ***EBEAM=BEAM ENERGY (MEV) > FORFAC = FORM FACTOR ( F( K ) ) 

C ***ETA=BEAM LENGTH/WAVELENGTH > BETA=VELOCITY/C( FOR MEDIUM) 

C ***GAMMA=RATIO OF BEAM ENERGY ELECTRON REST ENERGY 

C ***DIFFUN=D I FRACTION FUNCTION (I(U)) * DTHETA=THETA( DEGREES ) 

C ***THETA=THETA( RADIANS ) \ LENGTH=LENGTH FROM SOURCE ( METERS ) 

C ***RADFUN=RADIATION FUNCTION* MI DNU=BE GINNING FREQUENCY 

C ***CUE=BUNCH CHARGE IN COULOMBS MEREST = ELECTRON REST ENERGY 

C ***ERG=ARRAY FOR STORING VALUES OF RADIATED ENERGY 

C *** FOR CUMULATIVE ANGLES THETA 

C ***DNU= FREQUENCY INTERVAL *TERG=RADIATED ENERGY FOR A 

C *** SPECIFIC ANGLE THETA 

C 
C 

C ***ASSIGNMENT OF CONSTANTS*** 

CO=2 . 997925E8 
PI=3. 14159 
EREST = • 5117 
MU=1 . 2566E-6 
C 
C 

PRINT *, 'THE INDEX OF REFRACTION (AIR) USED IS 1.000268' 

PRINT *, ' DO YOU WANT A DIFFERENT VALUE?(Y/N)' 

READ (5,90 )ANSW 
90 FORMAT ( A1 ) 

IF ( ANSW .EQ.'N')THEN 
IND=1. 000268 
GO TO 100 
ELSE 

PRINT *, 'SELECT THE INDEX OF REFRACTION (REAL)' 

READ *,IND 
END IF 

100 CONTINUE 
C 

C DETERMINING THE SPEED OF LIGHT FOR THE MEDIUM 
C = CO/IND 
C 

PRINT *, 'ENTER BUNCH CHARGE (REAL) IN COULOMBS' 

READ *,CUE 
260 CONTINUE 
C 

CCC COMPUTATION OF CONSTANT (CALLED Q IN THESIS) 

C 



* 

* 



TRI00020 

TRI00030 

TRI00040 

TRI00050 
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CONST=MU*C*( CUE**2 . )/( 8 . *( PI**2 . ) ) 

CALL EXCMS ( ’CLRSCRN' ) 

C 

c 

PRINT *> ’ENTER THE DISTANCE FROM THE SOURCE (IN METERS) 1 
PRINT *, ’THE RECEIVER MILL BE PLACED* 

READ *> LENGTH 

PRINT *, ’SELECT THE FORM FACTOR DESIRED:’ 

PRINT * > ’ 1=UNITY > 2=TRAPEZOID 1 
READ *, I FORM 

IF ( IFORM .EQ. 1)THEN 
FORFAC=l . 0 
GO TO 290 
ELSE 

PRINT * y 1 ENTER THE TOP AND BASE VALUES (IN NANOSECONDS) FOR’ 
PRINT *,'FOR THE CURRENT TRAPEZOID FUNCTION' 

READ *>TOP>BOT 
C 

C THE NEXT 2 LINES CONVERT PULSE TIME TO DISTANCE 

BEE = C#TOP/2 . E9 
DEE = C*BOT/2 . E9 

END IF 

290 CONTINUE 

CALL EXCMS ( ’CLRSCRN’ ) 

C 

C BEAM ENERGY SELECTION: USED TO CALCULATE BETA A GAhMA 
C 

PRINT *, ’SELECT BEAM ENERGY:* 

PRINT * > * 1=50 MEV \ 2= SELECT YOUR OUT 
READ *,IBEAM 

IF ( IBEAM .EQ. 1 )THEN 
EBEAM=50 . 

GO TO 300 
ELSE 

PRINT *> ’INPUT THE DESIRED ELECTRON BEAM ENERGY IN MEV* 
PRINT *>’( INCLUDE DECIMAL POINT)* 

READ *>EBEAM 
END IF 

300 CONTINUE 

GAMMA=EBEAM/EREST 
BETA=SQRT ( 1 . - ( 1 . /GAMMA** 2 . ) ) 

CALL EXCMS ( 'CLRSCRN* ) 

C 

c 

PRINT *, ’SELECT A FREQUENCY BAND TO BE SOWED OVER AND THE’ 

PRINT *> ' INCREMENT IN MHZ: 1 = A SUM OF FREQUENCIES' 

PRINT *> ' FROM 10 MHZ TO 100 MHZ IN 10 MHZ INCREMENTS* 

PRINT *,'2 = SELECT YOUR OWN’ 

READ *» INCREM 

IF (INCREM .EQ. 1 )THEN 
BNU=10.E6 
D LOOP =9 
DNU=10.E6 

CALL EXCMS ( *CLRSCRN* ) 

GO TO 400 
ELSE 

CALL EXCMS ( ’CLRSCRN’ ) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
C CALL THE SUBROUTINE USED TO DEFINE FREQUENCY BAND OF INTEREST C 
C AND THE INCREMENTAL STEPS FOR SOWING C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL FREQ( BNU>DNU>DLOOP ) 

END IF 

400 CONTINUE 
C 

C OUTPUT GRAPH AXIS SELECTION 

PRINT *, ’SELECT OUTPUT TYPE* 

PRINT *,’1=ENERGY VS. FREQUENCY FOR A SELECTED ANGLE THETA’ 

PRINT *>’2=ENERGY( INTEGRATED OVER FREQ' 'S OF INTEREST) VS. THETA* 



TRI00100 

TRI00110 
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PRINT *, '3=F0RMFACT0R(F<K>) VS. FREQ. FOR A SELECTED ANGLE THETA' 

READ *, SELECT 

CALL EXCMS ( 'CLRSCRN ' ) 

IF (SELECT .EQ. 1 .OR. SELECT .EQ. 3) THEN 

GO TO 470 

ELSE 

GO TO 450 
END IF 

450 CONTINUE 

PRINT *, 'SELECT A RANGE FOR THETA AND THE INCREMENT AT WHICH' 

PRINT *, 'THE RADIATED ENERGY WILL BE CALCULATED' 

PRINT *, '1 = 0 TO 180 DEGREES IN ONE DEGREE INCREMENTS' 

PRINT *, '2 = SELECT YOUR OWN' 

READ *, NUMBER 

IF (NUMBER .EQ. 1 )THEN 
D2LOOP=181 
DTHETA=0 . 0 
DEGINC=1 
ELSE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C CALL THE SUBROUTINE USED TO DEFINE THE INCREMENTAL C 

C STEPS OF THETA C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL DEGREE ( DTHETA ,DEGINC ,D2LOOP ) 

END IF 

465 CONTINUE 

GO TO 500 

470 PRINT *, 'SELECT AN ANGLE (THETA) OF INTEREST.' 

READ *, DTHETA 

MIDNU = BNU 
TTERG =0.0 
D LOOP =D LOOP +1 
GO TO 510 
500 CONTINUE 

DO 600 K=l,D2LOOP TRI00140 

C 

C THIS OUTER LOOP WILL BE USED ONLY FOR OUTPUT TYPE 2 
TTERG=0 . 

MIDNU= BNU+( DNU/2. ) 

510 CONTINUE 

THETA=DTHETA*PI/180 . 

COSANG = COS( THETA) 

KONST = 2 . *SIN( THETA )/( ( 1 ./( IND*BETA ) ) -COSANG ) 

ATHETA = ( LENGTH*PI/C )*( ( l./( IND*BETA ) )-COSANG ) 

C 

C THIS INNER LOOP IS USED FOR ALL 3 OUTPUT TYPES BY INCREMENTING 
C FREQUENCY OVER THE CHOSEN RANGE 

DO 550 J=l,DLOOP 

IF ( IFORM .EQ. 1)THEN 
GO TO 515 

ELSE 

COSARG = 2 . *BEE*PI/C 
COSGAR = 2 . *DEE*PI/C 

KKONST = ( C**2 . )/( 2 . *( PI** 2 . )*( ( DEE** 2 . )-( BEE**2 . ) ) ) 
FORFAC=KKONST*( l./(MIDNU**2. ) )*( 1 ./( ABS( COSANG )**2 . ) )* 

D ( COS( COSARG*COSANG*MIDNU )-COS( COSGAR*COSANG*MIDNU ) ) 

END IF 

515 CONTINUE 

RADFUN=ABS( KONST*SIN( ATHETA»MIDNU )*FORFAC ) TRI00180 

TERG=CONST*(RADFUN**2. ) 

IF (SELECT .EQ. 1) THEN 

WRITE( 21,520 ) MIDNU, TERG 
520 FORMAT (2E 10. 5) 

ELSE IF (SELECT .EQ. 3 )THEN 

WRITE (21,525) MIDMJ,FORFAC 
525 FORMAT ( E10 .5,E10 .4 ) 

ELSE 

TTERG = TTERG ♦ TERG 

530 CONTINUE 
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END IF 

MIDNU= MIDNU+DNU 

550 CONTIMJE TRI00200 

IF( SELECT . EQ. 1 .OR. SELECT .EQ. 3) THEN 
GO TO 600 
ELSE 

GO TO 555 
END IF 

555 CONTINUE 

ERG! K )=TTERG* DNU 
GEE = ABS ( KONST ) ** 2 . 

WRITE! 21,560) DTHETA, ERG! K ) TRI00190 

560 FORMAT (2E10.5) 

DTHETA = DTHETA ♦ DEGINC 
600 CONTINUE 

C THE DATEND STATEMENT IS USED TO END THE DATA FILE 
C (NECESSARY IF USING DISSPLA GRAPHICS) 

DATEND=999 . 99 

WRITE! 21,660) DATEND , DATEND 
660 FORMAT (2F10.2) 

PRINT *, 'DO YOU WANT TO MAKE ANOTHER RUN?(Y/N)' 

READ ( 5,670 )RESPON 
670 FORMAT! Al) 

I F C RESPON. EQ . 'Y' ) THEN 
CALL EXCMS ( 'CLRSCRN ' ) 

GO TO 100 
ELSE 

GO TO 700 
END IF 

700 CONTINUE 

STOP TRI00210 

END TRI00220 

C 

c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 



SUBROUTINE FREQ! BFREQ, DNU, DLOOP ) 

REAL BFREQ, EFREQ,DIFF, DNU 
INTEGER DLOOP 

PRINT* ,' ENTER BEGINNING FREQUENCY IN MHZ! INCLUDE DECIMAL POINT)' 

READ* , BFREQ 

BFREQ=BFREQ*1.E6 

PRINT* ,' ENTER ENDING FREQUENCY IN MHZ! INCLUDE DECIMAL POINT)' 

READ* ,EFREQ 

EFREQ=EFREQ*1.E6 

PRINT* ENTER DESIRED INCREMENT IN MHZ! INCLUDE DECIMAL POINT)' 

READ* ,DNU 

DNU=DNU*1.E6 

DIFF=EFREQ-BFREQ 

DLOOP=INT ( DIFF/DNU ) 

RETURN 

END 



C 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

SUBROUTINE DEGREE! BDEG, DEGINC ,D2LOOP ) 

REAL BDEG, EDEG, DEGINC, DIFFF 

INTEGER D 2 LOOP 

CALL EXCMS ( 'CLRSCRN' ) 

PRINT *, 'ENTER STARTING AND ENDING ANGLES OF INTEREST! REAL ) ' 

READ *, BDEG, EDEG 

PRINT *, 'ENTER THE INCREMENTAL VALUE OF THETA! DEGREES ) DESIRED.' 

READ *, DEGINC 

DIFFF = EDEG-BDEG 

D2LOOP =INT(DIFFF/DEGINC)+1 

RETURN 

END 
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